Emerging understanding of the AI = 1/2 Rule from Lattice QCD 
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There has been much speculation as to the origin of the A J = 1/2 rule (ReAo/ReA2 — 22.5). 
We find that the two dominant contributions to the AI = 3/2, K — s> tttt correlation functions 
have opposite signs leading to a significant cancellation. This partial cancellation occurs in our 
computation of Re^ with physical quark masses and kinematics (where we reproduce the exper- 
imental value of A2) and also for heavier pions at threshold. For Re^lo, although we do not have 
results at physical kinematics, we do have results for pions at zero-momentum with m, — 420 Me V 
(ReA /ReA 2 = 9.1(2.1)) and ~ 330 MeV {ReA /ReA 2 = 12.0(1.7)). The contributions which 
partially cancel in Re^2 are also the largest ones in ReAo, but now they have the same sign and 
so enhance this amplitude. The emerging explanation of the AI =1/2 rule is a combination of the 
perturbative running to scales of 0(2 GeV), a relative suppression of Re^2 through the cancellation 
of the two dominant contributions and the corresponding enhancement of Reylo- QCD and EWP 
penguin operators make only very small contributions at such scales. 
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Introduction 

The 11 AI = 1/2 rule" remains one of the longest-standing 
puzzles in particle physics. It refers to the surprising fea- 
ture that in K — > tttt decays the final state is about 450 
times more likely to have total isospin 1=0 than 1=2. In 
terms of the (predominantly real) K — > tttt amplitudes 
Ao and A 2 , where the suffix denotes I, this corresponds 
to Re^4o/R e ^2 — 22.5. Perturbative running from the 
electroweak scale to about 1.5 -2 GeV contributes a fac- 
tor of approximately 2 to this ratio [TJ [5] ; the remaining 
factor of about 10 should come from non-perturbative 
QCD or, just possibly, from new physics. Lattice QCD 
provides the opportunity for the non-perturbative evalu- 
ation of A and A 2 , although it is only very recently that 
such direct K — > tttt calculations have become feasible. 
In this letter we summarise the emerging explanation of 
the AI = 1/2 rule from computations of ylo and A2 by 
the RBC-UKQCD collaboration. 

The first results from direct simulations of a kaon de- 
caying into two pions were presented in [3H5] . The deter- 
mination of Ao , where the two pions have vacuum quan- 
tum numbers, is particularly challenging and so far it 
has not been calculated with physical masses and mo- 
menta. We are striving to overcome technical issues such 
as the efficient evaluation of disconnected diagrams and 
the projection of the physical state through the use of G- 
parity boundary conditions [6 9 in order to evaluate Aq 
at physical kinematics in the near future. In the mean- 
time we have evaluated Aq and A2 for pions with masses 
of approximately 420 MeV [3] and 330 MeV [TO] at thresh- 



old, i.e. with the pions at rest. For these unphysical 
masses we do find a significant enhancement of the ratio 
ReAo/Re^2, albeit a smaller one than 22.5 (see the first 
two rows of Table [I]) . While investigating the origin of 
this enhancement we found a surprising cancellation in 
the evaluation of Re^2 , which significantly increases the 
ratio ReAo /ReA 2 . This suppression of Re^2 is the main 
result presented here. 

We have also evaluated A 2 with physical masses and 
momenta, obtaining a result for Re^ which agrees with 
the physical value and determining Im^2 for the first 
time [H [5] (see the third row of Table 1). In the eval- 
uation of Re^4.2 at physical kinematics there is a similar 
cancellation; indeed it is even more pronounced than at 
the unphysical masses in the first two rows of Tab. [I] 

In the next section we summarize the simulations we 
have performed, highlighting features of immediate rele- 
vance for the AI = 1/2 rule and referring to earlier pub- 
lications for other details. We then explain the partial 
cancellation of the two contributions to ReA 2 , which con- 
tradicts naive expectations from the factorization (vac- 
uum insertion) hypothesis. We also show that these two 
contributions have the same sign in Re^4o- We conclude 
by explaining how these features combine to provide an 
emerging understanding of the AI =1/2 rule. Of course 
a full quantitative explanation will require a calculation 
of ReAo at physical kinematics which is underway. 

Calculation of the Decay Amplitudes 

Our evidence is based on calculations from three Domain 
Wall Fermion (DWF) ensembles with 2+1 sea-quark 
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a' 1 [GeV] 


m„ [MeV] m K [MeV] ReA 2 [10" 8 GeV] ReA [10" 8 GeV] 


RcAq 
KcA 2 


notes 


16 3 Iwasaki 


1.73(3) 


422(7) 878(15) 


4.911(31) 


45(10) 


9.1(2.1) 


threshold calculation 


24 3 Iwasaki 


1.73(3) 


329(6) 662(11) 


2.668(14) 


32.1(4.6) 


12.0(1.7) 


threshold calculation 


IDSDR 


1.36(1) 


142.9(1.1) 511.3(3.9) 


1.38(5)(26) 






physical kinematics 


Experiment 




135-140 494-498 


1.479(4) 


33.2(2) 


22.45(6) 





TABLE I: Summary of simulation parameters and results obtained on three DWF ensembles. 



flavours (see Tab.[T]). Papers [H [3] describe a complete 
calculation of A 2 using the IDSDR (Iwasaki + Dislocation 
Suppressing Determinant Ratio) gauge action [TT] for (al- 
most) physical pion and kaon masses and realistic kine- 
matics. The ensemble was generated at a single lattice 
spacing a (a -1 ~ 1.4 GeV) chosen so that the volume 
is sufficiently large to accommodate the propagation of 
physical pions. In [3] a complete calculation of both Aq 
and A2 was carried out with the Iwasaki gauge action 
for 771 ~ 422MeV and m K ~ 737, 878 and 1117 MeV 
(here we present results for itlk — 878 MeV which cor- 
responds to almost energy-conserving decays) . Although 
the calculation was performed at threshold, this was the 
first time a signal for Re^4o bad been obtained in the 
direct evaluation of the K — > tttt matrix elements. A 
similar threshold calculation was presented in [lOj on a 
larger volume (24 3 ) with 777^= 329 MeV. The increased 
time extent of this lattice suppresses "around-the-world" 
effects in which one of the pions from the sink propa- 
gates in the forward time direction, crossing the periodic 
boundary and reaching the weak operator with the kaon. 
The calculation also used two-pion sources in which the 
single-pion wall sources are separated in time by a small 
number of time slices 6 (the results presented here are 
for (5 = 4). We find that this suppresses the (unphysical) 
vacuum contributions in the 1 = channel, significantly 
reducing the noise. In this way Re^o was resolved using 
only 138 configurations, compared to 800 in [3J. 

The amplitudes Aq and A 2 can be expressed in terms 
of the "master formula" 
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(7 = 0,2). (1) 



t = —V t * s V t d/V* s V u d and the Vij are elements of 
the Cabibbo-Kobayashi-Maskawa matrix. M 4 A/,lat = 
((tttt)j I Q\ &t I K) are the matrix elements calculated on 
the lattice. They are determined by fitting three-point 
correlation functions composed of a kaon source at t = 0, 
a two-pion sink at t = A, and one of the operators Q' at 
in the weak Hamiltonian inserted at all times < t < A. 
We fit the correlation functions Ci } i(A,t), 
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for <C t <C A, using a one parameter exponential fit 
to determine the matrix elements M J A ' f ' lat . All these cor- 
relation functions can be expressed in terms of the 48 
contractions enumerated in Section IV of [3] and labelled 
® through @. The contractions are functions of A and 
t, but we leave this dependence implicit, writing for ex- 
ample C 2)1 (A,t) =i v / 2/3{® + ©}._ 

The renormalization factors Z^^^ provide the con- 
nection between the bare lattice operators and those 
renormalized in the MS-NDR scheme at the scale /i, 

QY S (ri = Zlf->™(»,a)Q'}*(a). (3) 

The operators Qi on the left of ^ correspond to the 
conventional 10-operator "physical" basis, which is over- 
complete (see e.g. (H]). When calculating the renormal- 
ization factors, it is convenient to work in an equivalent 
"crural" basis of 7 linearly independent operators Qj with 
definite 5f (3)l x SU(3)r transformation properties (see 
eqs.(172)-(175) in [12]). z* (/j) + Ty» (fj,) are Wilson coeffi- 
cient functions. Fj is the Lellouch-Liischer factor relating 
the finite-volume Euclidean-space matrix element to the 
physical decay amplitude [T3] , 

Evaluation of ReA2= A2 receives contributions from 
the Electroweak Penguin (EWP) operators Qj and Qs as 



well as a single operator Q 
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where the superscript 3/2 denotes AI and the sub- 
script (27, 1) denotes how the operator transforms under 
SU(3)l x SU(3)r chiral symmetry. i,j are color labels 
and the spinor indices are contracted within each pair of 
parentheses. The subscript L denotes left, so that e.g. 
(s'tf = ~ 7 5 K) (^7^(1 - 7V). The 

AI = 3/2 components of the operators Qi,Q2,Q9 and 

3 /2 

Q10 are all proportional to Q,^ iy From all our simula- 
tions we confirm that the contribution from the EWP op- 
erators to Rc^2 is about 1%; e.g. for physical kinematics 
wc find RcA 2 = (1.381±0.046±0.258) 10" 8 GeV to which 
the EWP operators contribute -0.0171 10~ 8 GeV [1 [5] 
(the physical value is ReA 2 = 1.479(4)10~ 8 GeV). We 

therefore neglect the EWP operators in the following dis- 

3/2 

cussion. Chiral symmetry implies that 1) does not 
mix with the EWP operators so that ReA 2 is propor- 
tional to its lattice matrix element (the constant of pro- 
portionality is the product of the Wilson coefficient, the 
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Contraction ©. 



Contraction ©. 



FIG. 1: The two contractions contributing to Re^- 
They are distinguished by the color summation (i, j 
denote color), s denotes the strange quark and L that 
the currents are left-handed. 
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FIG. 2: Contractions ©, -© and © + © as functions of 
t from the simulation at physical kinematics and with 
A = 24 . 



renormalization constant, finite- volume effects and kine- 
matical factors; see [5] for a detailed discussion). 

Fierz transformations allow the K — > tttt correlation 

3/2 

function of Q/ 2 >j u to be reduced to the sum of the two 
contractions illustrated in Fig.[T] labeled by ® and @. 
The two contractions are identical except for the way 
that the color indices are summed. A 2 is proportional to 
the matrix element extracted from the sum ® + © . The 
main message of this letter is our observation from all 
three simulations that ® and © have opposite signs and 
are comparable in size. This is illustrated in Fig.[2]for the 
results at physical kinematics from [H [5] , where we plot 
©, -© and © + © as functions of t. We extract A 2 by 
fitting © + @ in the interval t € [5, 19] where there is a 
significant cancellation between the two terms. A similar, 
although not quite so pronounced cancellation occurs at 
threshold for physical masses and for the heavier masses 
studied in [31 [TO], see Fig. [3] for example. 

It has been argued that the factorisation hypothe- 
sis |14j works reasonably well in reproducing the experi- 
mental value of Ai (see e.g. Sec.VIII-4 in [T5]). In this 
approach, the gluonic interactions between the quarks 
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FIG. 3: Contractions ©, -© and © + © as functions of 
t from the simulation at threshold with ~ 330 MeV 
and A = 20. 



combining into different pions are neglected and A2 is 
related to the decay constant f n and the form factor 
close to zero momentum transfer. On the basis of color 
counting, one might therefore expect that © ~ 1/3®, 
whereas, for physical kinematics, we find © ~ —0.7® 
and that nevertheless ® + © leads to the correct result 
for A2. Thus the expectation based on the factorisation 
hypothesis proves to be unreliable here. 

Following the discovery that ® and © have opposite 
signs we examined separately the two contributions to 
the matrix element {K°\(sd)[ J (sd)L\K a ) which contains 
the non-perturbative QCD effects in neutral kaon mix- 
ing |llj . The two contributions correspond to Wick con- 
tractions in which the two quark fields in the K° interpo- 
lating operator are contracted i) with fields from the same 
current in (sd)L(sd)L and ii) with one field from each of 
the two currents. Color counting and the vacuum inser- 
tion hypothesis suggest that the two contributions come 
in the ratio 1:1/3, whereas we find that in QCD they 
have the opposite sign. 

We stress that it is only the correlation function ®+© 
which has a time behaviour corresponding to £'( 7r7r ) 2 . 
Because the calculation is performed in a finite-volume 
-^(7T7r)2 7^ ^(tttt)o an d © and © individually have an 
isospinO component. If E/ nn \ 2 = ttik then ® + © is in- 
dependent of t away from the kaon and two-pion sources, 
and this is what we observe, particularly in Fig. [2] where 
the energies are matched most precisely. 

We postpone a discussion of the implications of these 
results to the AI =1/2 rule until the next section, but 
we believe that the partial cancellation observed in the 
evaluation of Ai is a significant component. 

Evaluation of ReA : The evaluation of A at phys- 
ical kinematics has not yet been completed. The results 
presented here are obtained at threshold, with the two pi- 
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ons in their zero-momentum ground state with each pion 
at rest up to finite-volume effects. Even at threshold 
we have had to overcome many theoretical and technical 
problems, including the evaluation of the 48 contractions 
contributing to the correlation functions, the renormal- 
ization of the operators in the effective Hamiltonian, the 
subtraction of power divergences and the evaluation of 
the finite-volume corrections. The threshold calculations 
do not require however, the isolation of an excited state. 
The pions in a physical decay each have a non-zero mo- 
mentum in the center-of-mass frame, which corresponds 
to an excited state in lattice calculations. Given the poor 
statistical signals after the subtraction of power diver- 
gences and the evaluation of disconnected diagrams, the 
evaluation of A Q at physical kinematics is currently im- 
practicable with standard techniques and is the main mo- 
tivation for our development of G-parity boundary con- 
ditions [SHU]. 

With the two pions at threshold we find [31 [TO] 

ReA _ ( 9.1(2.1) for m K = 878 MeV, m, = 422 MeV 
ReA 2 ~ 1 12.0(1.7) for m K = 662 MeV, = 329 MeV. 

(5) 

While these results differ significantly from the observed 
value of 22.5, because the calculations are not performed 
at physical kinematics, there is nevertheless already a sig- 
nificant enhancement in the ratio and it is interesting to 
understand its origin. In Tab. [Til we present the contri- 
butions to KcAq from each of the lattice operators in the 
24 simulation with a -1 = 1.73(3) GeV and from each 
MS-NDR operator at renormalization scale 2.15 GeV. In 
both cases, the dominant contribution comes from the 
current-current operators Q 2 - 

Since in a finite- volume E^^ 2 ^ E(- w \ , one cannot 
satisfy the condition ttlk — E n7r for both isospin chan- 
nels simultaneously with the same quark masses. Here 
we quote results using the fixed meson masses quoted 
in Eq. ([5]), which is sufficient for our current discussion. 
For these masses E ( ^ )o = 766(29) MeV (629(15) MeV), 
E { ^ )2 = 876(15) MeV (668(11) MeV) for the 16 3 (24 3 ) 
lattice. A study that interpolates in the kaon mass to 
make both decays energy-conserving may be found in [3] . 

The dominant contribution from the lattice operator 
Q 2 to the AI =1/2 correlation function is proportional 
to the contractions 2-® — @ and corresponds to typel 
diagrams in the language of [3j (see Fig. 3 in [3j). In Fig. [4] 
we show the total contribution of Q 2 to the correlation 
function, as well as the total connected contribution and 
that of typel diagrams given by ^|{2-@ — ©}. The 
errors on the total contribution are dominated by the 
disconnected diagrams. The observation that © and © 
have opposite signs leads to an enhancement between the 
two terms rather than the suppression in the factorization 
approximation © = Similarly, in the case of Q\, 

the typel combination ^{2-© — ©} is dominant. In 
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TABLE II: Contributions from each operator to ReA$ 
for m K = 662 MeV and m T = 329 MeV. The second 
column contains the contributions from the 7 linearly 
independent lattice operators with 1/a = 1.73(3) GeV 
and the third column those in the 10-operator basis in 
the MS-NDR scheme at fi = 2.15 GeV. Numbers in 
parentheses represent the statistical errors. 
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FIG. 4: Contributions of Q 2 at to Re^4o (purple crosses). 
The blue squares and black circles denote the connected 
and typel contractions respectively. 

this case both the correlation function and the Wilson 
coefficient Zi(fi) +ryi(fJ,) are negative, so that the overall 
contribution adds to that from the correlation function 
of Q 2 . 

Finally we note that in our data ReA 2 shows a much 
stronger mass dependence than ReAo, which was also 
expected in SU(2) chiral perturbation theory [Jj)j. We 
attribute this to the partial cancellation between © and 
© in Re^2- Our results for ReA 2 and ReA are given in 
TabHJ 

Conclusions 

From our recent computations of K —¥ rnr decay am- 
plitudes a likely explanation of the AI = 1/2 rule is 
emerging. In particular, we find that in the evaluation of 
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Re^2, which is proportional to the sum of two contrac- 
tions ® + @, there is a significant cancellation between 
the two terms. The naive expectation based on the fac- 
torization hypothesis suggests that © « |®, whereas in 
QCD we find that they have the opposite sign. (The two 
terms contributing to Bk similarly have opposite signs 
contradicting expectations from the vacuum insertion ap- 
proximation.) 

The evaluation of A$ at physical kinematics has not 
yet been performed. Our simulations at threshold with 
772 = 329 MeV and 422 MeV show that the dominant 
contributions to A$ comes from the current-current op- 
erators, with only small corrections from the penguin op- 
erators. This is true whether we express the results in 
terms of the bare lattice operators at a -1 = 1.73 GeV 
or the MS-NDR renormalized operators at /i = 2.15 GeV 
(see Tab.[i"H). Although 48 contractions contribute to the 
1 = correlation function, in our simulations the largest 
contributions again come from contractions ® and @ 
with relative signs which enhance Re^4rj- 

References to estimates of the amplitudes using ana- 
lytic or model approximations are presented in the re- 
views [TTJ [TH] . We note that a suppression of Re^2 and 
an enhancement of KbAq was found in [19] using the 1 /N 
expansion with a particular ansatz for matching the short 
and long-distance factors at scales 0.6-0.8 GeV. 

The results presented above indicate that Re^2 is very 
sensitive to the choice of quark masses and momenta; a 
sensitivity we attribute to the partial cancellation of the 
two contributing contractions. On the other hand, there 
is no such cancellation in ReA and indeed the results 
depend much less on the masses and the values we find 
are already close to the experimental result. Of course 
before we can claim to understand the AI =1/2 rule 
quantitatively, we need to reproduce ReAo/ReA2=22.5 
at physical quark masses and kinematics and we are 
currently undertaking this challenge. Nevertheless, from 
the results and discussion of this paper it appears that, 
in addition to the well known perturbative enhancement 
of ReAo/Rev4.2, the explanation is a combination of 
a significant relative suppression of Re^ as well as 
some enhancement of ReA with penguin operators 
contributing very little. 
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